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Abstract 

Two new general representations (the series and the integral) for the mass current 
j in weakly inhomogeneous superfluid A-phase of Helium-3 are obtained near zero 
of temperature by solving the Dyson-Gorkov equation. These representations result 
in additional correcting contribution to the standard leading expression for j which 
is of first order in gradients of the orbital angular momentum vector I. The total 
supplementary term is found as integral, and, provided the London limit holds, 
the procedure is advanced to expand it at T = asymptotically by the Laplace 
method in powers of gradients of I. Three special static orientations of I with 
respect to its curl are considered to calculate the higher correcting terms up to 
third order. Coefficients at the quadratic terms are estimated numerically, new 
cubic contributions are found which contain the logarithm of the London parameter. 
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1 INTRODUCTION 



Superfluidity of Helium-3 is firmly in focus of intensive theoretical and experimental stud- 
ies [1]. Although considerable attention is paid last years to such problems as quantized 
vorticity and interfaces [1-4], the weakly inhomogeneous A-phase of Helium-3 ( 3 He — A) 
also can be a subject for theoretical investigation. This phase originates due to p-wave 
spin triplet -BCS-pairing [5, 6], and it demonstrates rather unconventional behaviour [7]. 
Its peculiarity can be seen, for instance, from the mass current j which is acknowledged 
to be of first order in gradients [5-7] : 

Jo = pv s + ^rot {pi) + j an , (T = 0) (1) 

Jan = ~ 7T~ C iff ■ rot 

where the first two terms are habitual for a nodes-free p-wave superfluid, while the famous 
anomalous j an witnesses the unusual manifestations of the ground state of 3 He — A. The 
point is that j an is caused by nodes existence in the gap on the Fermi surface for real 
3 He — A ]7, 2]. In (1) p is the liquid density, m is the atom mass, v s is the superfluid 
velocity, I is the weakly inhomogeneous orbital angular momentum vector (hat implies 
a unit vector) and Co ~ p. The Eq.(l) has been deduced by many authors in different 
approaches: by solving the Gorkov [8-10] or the matrix kinetic equations [11], as well as 
in [12] directly with the use of the ground state wave function. 

Alhough the Eq.(l) and the corresponding physical picture have been broadly dis- 
cussed [13-21], slight indications can be seen in [8, 12, 20] that higher corrections to (1) 
might occur thus causing difficulties of the superfluid hydrodynamics at T = 0. Indeed, in 
the course of a phenomenological consideration of the free energy of 3 He — A Volovik and 
Mineev have found one of these corrections in the form Xorb D l a d l a , where D — dt + v s • d 
[13]. Aiming to check on their own the presence of j an in (1), Combescot and Dombre 
have developed a microscopic calculation [10] which has allowed to claim at T = the 
quadratic correction \l x rot I |(rot l)± in the current perpendicular to I, and the terms 
| fx rotf| (v s — (l/4ra)rotZ^ and | fx rot f| (9^2 + d 2 fi) in the current parallel to T. 

Alhough [10] has eventually been aimed at the case T = 0, and the Gorkov equation 
has been solved exactly after linearization of the order parameter, an intermediate high- 
temperature approximation to obtain a manageable formula for the ^-integrated Green 
function has not been avoided. The difference between the approximate and exact formu- 
las has been considered as responsible for the second order corrections at T = 0. However, 
due to the strategy adopted in [10], any chance to be accurate with T — > and to ben- 
efit from the exact solution of the governing equation has been lost. Consequently, the 
coefficients at the quadratic corrections and some other formulas have become artificially 
complicated, e.g. requiring additional efforts to state that they are finite and nonzero. 
As to the higher contributions at T = 0, the regular expansion procedure has not been 
convincingly formulated (see discussion in [23]). 

Since unambiguous procedure of asymptotic expansion of j which would lead to a more 
deep knowledge of the ground state of 3 He— A seems a meaningful technical problem, a way 
to overcome mathematical difficulties of [10] has been advanced in [22, 23]. That is another 
method to solve the Dyson-Gorkov equation has been chosen which gave immediately 

— * 

new representations for the fermionic Green functions and thus for j. It is well known in 
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mathematical physics that Green function of a Sturm-Liouville operator can be written 
either as integral or as series in eigen-functions. The second possibility is meant here. 
After the subsequent paper [24] it became clear that this way looks properly because the 
new representations for j admit T — > accurately and can rigorously be studied by the 
Laplace method (steapest descent) provided the characteristic length in the texture is 
much longer then the coherence lenght £o (the London limit). As the result, it is possible 
to deduce corrections to (1) systematically in the form of asymptotic series in powers of 
gradients of I. 

This paper completes [22-24] and is organized as follows. Sec. 2 contains the outline of 
the problem which is the same as in [10] (i.e. the approximations and notations). Sec. 3 is 
concerned with the solution of the ordinary non-homogeneous differential equation related 
to the Dyson-Gorkov equation and with the calculation of the mass current in the form of 
series. Sec. 4 contains the integral representations for that series as well as various limits 

— * 

for the correspondingly written j: lowest "gradient" limit and zero temperature limit. 
Three particular cases of mutual orientation of rot / and / are considered in the Sec. 5 to 
obtain corrections to jo (1) & t T = up to third order. Apart from the quadratic terms 
predicted in [10], new cubic contributions are found which contain the logarithm of the 
London parameter. The numerical coefficients at the second order terms are provided. 
Discussion in the Sec. 6 concludes the paper. The present investigation supports rigorously 
the early results of [10, 13] and would be useful for any systematical microscopic approach 
to correct observables in 3 He — A. 



2 THE PROBLEM 

Since our main task is to calculate the mass current j by means of normal Green function, 
let us start with the standard matrix Dyson-Gorkov equation: 

d T g(k, k') - j d 3 k"H(k, k") g(k", k') = (2nf5 {3 \k - k') 5(t - r'). (2) 
Here r is "imaginary" time, g{k, k') is the 2x2 matrix of normal and anomalous two-point 

— * — * 

Green functions and H{k, k") has the form: 

H(k P>\ = ( ^" 5(3) ( k - k ") (2^)- 3 A(^, k") \ 
l ' ] \(2n)- 5 A*(k",k) -foS^^-h"))' 

where = (k 2 — hp) /2m, &f is the Fermi momentum and A(k, k") is the order parameter 
of 3 He — A. We shall calculate j by the formula 

i^- 1 ^)- 3 /^^^. (3) 

As far as the Refs. [22-24] have been conceived as technical improvement of the Ref. [10], 
the framework (i.e. approximations, notations) turns out to be unaffected here, and one 
should be referred to [10] for certain details. 

It is appropriate to re- write (2), (3) in the mixed coordinate-momentum representation 
[12, 25, 26]: 

H{k,r) = (2tt)- 3 / d 3 qH{k + q/2, k-q/2)e i ^, 
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g^r) = (2ny 3 J d 3 q g(k + q,k) e'** , (4) 

where r — \(fi + r*2) is the center of mass coordinate and the momentum k is conjugated 
to fi — f 2 . For instance, it is not difficult to check that the relation 

(2tt)- 3 J d 3 k" A(k, k") g(k", k') = J d 3 r e * A(k - iB r - % -d y , y) g % ,(r) 

holds by (4). Applying (4) to (2) we set in the lowest in 5/E F order: 

€(k-id) ~ 6 -im^k-d, 



y= r 



(-> -> i -> \ 
k-id r - -d y ,yj 



« A(k, r) = 5(k ■ Ai(f) + ik ■ A 2 (r)) , (5) 

y=f 



where 5 is the gap amplitude, k is unit reciprocal vector and the orbital momentum 
vector is given by A x x A 2 = /. The Eqs.(5) read that (2) can be written in the mixed 
representation in the form: 

f£-ic F k-d A(k, r) \ 

where g = g^(f), £ c F (k — k F ), c F is the Fermi velocity, and uj is fermionic Matsubara 
frequency. 

Thus we have obtained the approximate Eq.(6) which can be nicely treated as 1- 
dimensional because the spatial differentiation is along the directions labeled by k. Indeed, 
in [26] a gradient expansion method is presented to study dynamics of spatially inhomoge- 
neous systems provided inhomogeneities are slow compared to the relevant length scales. 
As the result, a separation of 3-dimensional problem into a collection of 1-dimensional 
subsystems occurs. Proofs useful for justification of our approach can be picked up from 
[26]. ^ " " ^ ' ^ 

As far as we are interested in j in arbitrary point, say, O, let us define the spherical 
coordinates p, 9, centered at it and linearize the slowly varying order parameter as 
follows: 

A(k, r) « A(k, p = 0) + ap = a(p + p) + iA , (7) 
A = lmA(k,p = 0) . 

As the physical result in O is assumed to be independent on the choice of the point, it 
can be calculated at any r with r — > O in final formulas. Therefore, we shall solve (6) at 

^ — * 

f = pk so that k ■ d is simply d/dp and put p = in the result [10]. 

To be precise, we shall consider our problem for the coherence length £ much smaller 
than a length of the orbital vector / variation: 

£o = y « \d®T\-\ 

or 

\ = £o\d®T\ < 1 (8) 
X 

(the London limit). The parameter a (7) depends on the angle variables, vector's v s 
components and first derivatives of / taken in O, and it will be written explicitly in the 
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Sec. 5. From (8) it is seen that the condition ap < 5 ensuring the linearization (7) implies 

p/£o Z X 2 i and holds better provided x 2 is greater. 

Changing the variable x = (a/c F ) 1 / 2 (p + po) and eliminating £ from the L.H.S. in (6) 
one gets: 

(iu + H)G = e^-F)- 172 !, (9) 

where 

H = iJac F a 3 — Jac F a x x + Aa 2 (10) 

dx 

is the Hamiltonian and <7j are the Pauli matrices. In this case (3) becomes 

;= ^(stt 3 ^)- 1 ! an% (p-^jj , (n) 

where J is the ^-integrated Green function: 

J(x) = | d^e-^^^ 172 G n (x), (12) 
and G(x) is to be determined from (9). 

3 SOLUTION OF THE NON-HOMOGENEOUS EQ- 
UATION IN FORM OF SERIES 

To solve (9) let us take G(x) in the form 

g ^(';: /:)■ ->)• (13) 

where /ii i2 = ^i,2( x ) and /i )2 = /i,2( x ) are now to be determined. Adjoint action 
u _1 <7iU = a 2 (cycl.perm.) of the unitary matrix u on the Pauli matrices transforms 
H (10) to H em : 

u^Hu = H em , n em ^(_.^_ a+ < ^°-) > (14) 

where a ± = x =p d/dx. The operator Ti, em reminds the Hamiltonian of spinning electron 
in constant homogeneous magnetic field. It is straightforward to obtain its eigen-values 
E , ±E n and eigen-f unctions § , $t ( n > 1) I 22 ]'- 

$ °=U(x))' £ » = - A ' 

Us) ( VEn + sA Vn-l(x) 



" I -«A - sA Vn(x) J ' SE "' (15) 

where s = ±, E n — ^ A 2 + 2ac F n and ?/> n (x) are the Hermite functions. 
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Let us use (14) to pass from (9) to the equation 

h 



fj -^-x')e^)- 1/2 (^ ; , 



(16) 



to calculate Gu. The Dirac 5-function is placed in the R.H.S. of (16), the unknown h, f 
depend now on x, x' and the required entry is given by 



G n (x) = /dx / (/ i (x,x / ) + /(x,x')) 



(17) 



To solve (K3) i1 is natural to expand ( ^. ) : the functions (lo): 



(x,x') = 5(x')§o(x) + £ E^M ■ 



(18) 



s=± n=l 



We calculate 5(x'), &i s) (x') using orthogonality of the vectors (15) [22, 24], and from 
(12), (17), (18) obtain J: 



J = 7i^ac f 



<£q,^0> + y ^ 

ioj + E 4l r< iu + sE n 



(19) 



where < •,• > stands for Hermitian scalar product. The representation (19) for the £- 
integrated Green function is alternative to that which has been found in [10] (the Eq.(34)) 
as a quadratic combination of parabolic cylinder functions. 
Now summation over uj is straightforward [27] and one gets: 



ac z 



1 oo 

n( Wo 2 + ? E (>n-i + ^n) + 



n=l 



/ n=l ^ n 



(20) 



where n(e) is the Fermi weight. Inserting (20) to j (11) we obtain the required general 
representation for the mass current near zero temperature [23, 24]. Due to the explicit 
dependence on /3, the Eq.(20) admits T — > as well: one has to replace n(E ) by the 
Heavyside function 9(E ) and tanh(/5£'„/2) by 1. It can be argued that the second term 
in (20) is irrelevant with regard to the angle integration and therefore j acquires the final 
form as the series: 



/r a 00 

dnky/ac; n(E )^ + - E (^'-i - € 

L A n=l 



tanh(/3£ n /2) 
E„ 



(21) 



4 INTEGRAL REPRESENTATIONS AND THEIR 
LIMITING CASES 



4.1 INTEGRAL REPRESENTATIONS 

Practically, it is more convenient to be concerned with an integral representation equiva- 
lent to (21). Such representation has been found in [23] at zero temperature, and here we 
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shall deduce it for (3 1 Y.J and j at T ^ [24]. These representations will admit special 

limiting cases. 

First of all, one has to rearrange (19) [22]: 



oo 



v -, 9 / iu — A iu + A \ . . 

By the formula 

/T^M^ + e 2 )- 1 = -1/2 

u> 

we see that the odd in oo part of (22) is responsible for the second term in (20) and it is 
enough to consider only the "even" part of (22): 

A OO 

Je = x -== E < [W + n + I)" 1 - (|A| 2 + n)^} , (23) 

Z V " C F n=0 

where |A| 2 = (c^ 2 + A 2 )(2ac F )- 1 . 

By the formula (AI.2) (APPENDIX I) the series (23) can be expressed as the integral 



/ \ 1/2 oo 

J e = -A J dt^anht) 1 /^- 2 ^^- 2 ^! 2 * 



(24) 



Thus one can go further: 

1/2 



(\ l/Z oo 
— j /dt(tanht) 
OLC F J J 



1/2 





x (T0 2 (O, ir)) e -x 2 tanh^(AVa CF )t ? (25) 

where the elliptic theta function i? 2 [28] implies the series 



oc 



E« (m+I)2 = ^(0^r), 

m=0 Z 

r = (— l/7r)loga, and a = exp(— 4n 2 T 2 t/ac F ). Changing the integration variable t i — > 
Kt, k — ac F ((3/2) 2 , one can rewrite (25) more suitably for studying the limiting cases 
below: 

1/2 





„2 , 



\ t ) 

x e(t)e- x2tanh(Kt) - (A/3/2)2 ', (26) 

where G(t) = (Tit) 1 / 2 tf 2 (0, i nt). 

The Eqs. (25) and (26) are just to be substituted to (11) to get the general integral 
representations for j near zero temperature. These representations are very convenient in 
calculating higher corrections to (1). Before to proceed to it in the Sec. 5, let us consider 
some particular limits for (25), (26). 
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4.2 LIMITING CASES 

Let us represent the situation by the following "commutative diagram" : 

" — 1 T 

(27) 

II -i+ (1) , 

where the horizontal arrows 1,3 mean T — > and the vertical ones 2,4 - the lowest 
"gradient" approximation. 

To begin with, the usage of the limit 

lim tf 2 (0,rr) = r~ 1/2 
[28] in (25) allows j to be written at T = as follows (the arrow 1): 

j = -3p (87rc F ) _1 JdttkAF (x 2 , A 2 /a c F ) , (28) 



where p = k^/3n 2 (two spin projections are taken into account), and the function F(p, q) 
is given by 

F(p, q) = Jdt f^l\ e -Pt^t- q t ^ q > 
o ^ ' 

To perform the lowest "gradient" approximation (the arrow 2), we replace tanh t by 
t in F(p, q) (see (8) and (AII.6)) and (28) takes the form: 



j = -3p(87rc F )- 1 | dQkA | dte _ V x2+ ^Fj 



o 

aA 



= -3p(8n)-' fdQk - (29) 

J ac F x z + A^ 

(compare with the Eq.(44) in [10]). To understand (29), let us recall the representation 
for j which has been discussed in [19, 20]: 

j = 3p (87T)- 1 J dVl k (^(k ■ d) arctan (^^j ) , ( 30 ) 

and which results in (1) with I = A x x A 2 and (v s )i = 2~ l A x ■ <9;A 2 . Using 

"a" d £;-A 2 A 



k-d 



c F dx ' fc • Ai x v /ac; ' 



it is easy to check coincidence of (29) and (30). So without any special gauge for a it is 
seen that (1) is the lowest London limit approximation to (28). 
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To consider the steps 4 — > 3, one should replace tanh(/tt) by Kt in (26) (the arrow 4) 
due to steepest descent validity at A ^ for large (5 : 

p-'^Je = -^P 2 J dte(t)e-« A WV 2t . (31) 

w 

The R.H.S. of (31) is the Laplace transform of and can be expressed through the 
so-called Yosida function Y [29]: 

oo oo 

a 2 J dt 9(f) e- 2 * = 1 - / cogh2 ^= = 1 - Y(a). (32) 

To check (32) it is enough to integrate its L.H.S. as the series and to re-express the answer 
by the Poisson summation formula, whereas cosh _2 (y) in the R.H.S. has to be expanded 
first in exp (— 2y) and then integrated [24]. Using (32) one obtains: 

J=-3p(8,)-7d fi ^(l-Y(^)) (33) 

(the point 11 on (27)). The Eq.(33) is just the leading j "dressed" by thermal corrections 
which has been found by Cross [8]. At T = (the arrow 3) Y(oo) = and we recover (1). 



5 EXPLICIT CALCULATIONS 

This section is devoted to the main problem of the present paper. That is it will be 
concerned with asymptotical expansion of j (28) in order to deduce the London limit 
corrections to (1). At fixed k the overall phase of the order parameter A(k, r) can always 
be changed to make a (7) a real positive. Thus the Eq.(7) can be thought of as 

exp(-iip) A(k, r) = A + ap, (34) 

where 

a = 5M exp(i(7r/2 - ip)) . (35) 
A = 5 sin 8 exp(i(0 — ip)) . 
In (34), (35) the phase ip is to be adjusted while M is given by 

M — - cos 2 6 d 3 T 2 + 2m sin 2 9 e 1 ^ (vi cos + v 2 sin 0) + 

i ~ 1 - 
+ sin 9 cos 9 e l<p (2mv 3 + — div I - — I ■ rot I) + 

+ 1 sin 9 cos 9 e-* {-d^ - d 2 h + i&h ~ d 2 T 2 )) , (36) 

where v = v s (APPENDIX IT). Without loss of generality A 2 (0) can be chosen along 
I x rot I so that d%l\ , become zero and thus div I = d\l\ + d 2 l 2 [10]. Besides, / • rot I is 
d\l 2 — d 2 l\ once the third axis is chosen along 1(0). Moreover, dil 1: d 2 l 2 can be excluded 
from the consideration [10]. Therefore, apart of v\ and v 2 , only d\l 2 + d 2 l±, 2mv 3 — 
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(1/2)/ • rot/, and ^3/2 = rot/ x / are the relevant gradient combinations. Besides, no 
difference is expected once M. is considered as dependent separately either on <9i/ 2 + d 2 h 
or 2mt> 3 — (1/2) /-rot/. 

To make calculations manageable it is appropriate to put a part of gradients in M. 

— * 

equal to zero so to consider the dependence of j on the remaining ones. Clearly, it is 
not necessary to enumerate all the possible cases, but it is enough to point out rather 
characteristic combinations. To this end, let us take M. in the following reduced form: 

M — — d 3 T 2 cos 2 9 + (2mv - dfy sin 9 cos 9 

' cos 2 6 + -^sin# cosfle^ (37) 



(at d{l 2 + d 2 h — 0, v = v 3 ). Once ip is obtained explicitly so that a G M + we get: 

5 2 sin 2 9 2 / sin \ 2 3 
aCp= Q ' X ° = U 2 tan 2 #J g ' 



6 Q= cos(, U!ta?-^ • (38) 



where 



1 / ! sin n \ /I cos o i i 



Q 2 \x! tan6y \xi tan# %i tan 2 9) 
A convenience is apparent after [24] to integrate by parts in F(p, q) in (28) so that 



Q 

with 



F(x 2 , A 2 /ac F ) = 1 (l + $(x 2 , Q)) (39) 



00 



$(X 2 , Q) = / e-O* ( ./^! e x 2 (t-tanh t ) I ^ (4Q) 







We have used in (39), (40) the following circumstance. According to the Eq.(34), the 
order parameter in O is A and the square of its modulus has the simple form: 

|A | 2 = |A| 2 = A 2 + ac F x 2 = 5 2 sin 2 9 

by the Eqs.(7) and (35). Now, from (38) it is seen that Q in (39), (40) is just \A\ 2 /ac F . 
The first term in (39) is responsible for j (see (29)), while the second one - for the total 
correcting contribution: 

Jcorr = - 3 P (87T)- 1 f dfi k $ ( X 2 , Q) . (41) 

In what follows we shall investigate (41) with M. (37): 

Jcorr = 3p(87re )- 1 / dtl k cos 9 (j^ - ^j^) $ (x 2 = 4 , Q) ■ (42) 

Varying xi, X2 in (42) the following three cases (Examples 1,2, and 3) can be obtained. 
Fixing X2 (or Xi) ^ 1 an d tending Xi ( or Xi) to infinity we shall get Example 1 (or 
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Example 2). Taking Xi = X2 = X ^ 1 we shall arrive to Example 3. Each time our 

— # 

attention will be called to the quadratic and cubic contributions to j cor r, i-e. to the terms 
proportional to (^oX 4 ) -1 an d (■CoX 6 ) 1 - Clearly, the case 1 corresponds to rot / parallel to 
/ and the case 2 — to rot I perpendicular to /. Therefore, the case 1 implies all the three 
contributions in (1), while the second one corresponds only to the pure orbital content of 
(!)■ 

5.1 EXAMPLE 1: rot I IS PARALLEL TO I. 

Here only 2mv — d\l2 7^ in M. (37). As far as we deduce from (38) that 

Q = x 2 1 tan #l , x o = , 

where (CoX 2 )' 1 = 2mt) — di/2 > 0, the Eq.(42) reads only the third component, say, j to 
be nonzero now: 

00 

3p 1 f 2 udu 







where 



and u = I tan 9 \ . Concrete form of T is not of importance for us. It is enough to know 
that T{s) — > const as s — > 0, and 

^u) = ^+^ + ..., x \»l, (44) 

where a= —1/3. 

To do the estimation let us break the integral over u into two parts: 

j/^(^TW = Ul + u ^ < 45 » 

where 



00 00 
f _ 9 . wd-u a /" . o r/n dw 



because x 2-u ^> 1 is valid, and we are interested in contributions of total degree in x n °t 
less than —6. As to Lq, 

x 2 

1 f ^ wdw 



'' = ?/^(T 





x 4 / (1 + «7x 4 ) 5/2 ' 

its denominator can be expanded in powers of w 2 /x 4 so that U\ will acquire the form of 
series where each term is given by the appropriate integral. The upper bound x 2 °f 
these integrals can be extended to infinity provided the integral is convergent, either some 
regularization by counter-term is needed. It is easy to see from (44) that / Tu du means 
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only the logarithmic divergency, / Fv?&u — both the quadratic and logarithmic ones and 
so on. 

Let us do the first subtraction writing Lq — X + Y , where 



V'l HU) + / - 1 ) MdM ' Y = ^jr(u)udu. (46) 



X 



V / \ \ V / / V 

o \\ *■ / ) * o 



Once the integral Y is divergent logarithmically at x 2 ~~ 00 ) it can be represented ap- 
proximately: 

1 °° 

X x o u + 1 / 

Let us turn to X (46). The asymptotic (44) tells us that X is divergent at x ~^ 00 > 
and the whole (44) is needed for regularization. It is not difficult to realize that the total 
contribution of the order x~ A appears as that counter-term which results from X once T 
is replaced by a/w 2 . Therefore, 

C/ 2 +X = ^ (-| + log 2) . 

Therefore at a = —1/3 

" + * = sj? (! + / ( 3 * * M) + iTTl) d " + log 2?) • (47) 

Finally, the use of the Eqs.(43), (45) and (47) enables the third component of j to be 
completly written as follows: 

*=T6M i+ ? log S' (48) 



where 



4 °° 1 

i„ g2e = - + /(3^ M + — )d«. 





All the quadratic corrections predicted in [10] are zero in the present gauge, and the lowest 
one turns out to be cubic with the logarithm of the London parameter. 

5.2 EXAMPLE 2: rot I IS PERPENDICULAR TO I 

In this case the Eqs.(38) result in 

Q = (xtan^) 2 , x = — x| tan^| sin , 

where (£oX 2 ) -1 = d^h > . From (42) the components 2 and 3 of j cor r are zero, whereas 
the first one acquires the form: 

3p 1 f fT , l 2 2\ udu 



f o 9x udu , . 

4 ^I^WTW' (49) 
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where 

'tanht „ „ 1 



o 

■u = | tan#| and the formula 

i 



J dv y/l=tfe** = j 1 F 1 ( i 2; p) , (50) 
o 

p = M 2 (t- tanht) > 

is used to re-express the integration over v = sin0. In (50) \Fi is the Kummer func- 
tion [28]. The relevant analytical properties of T are the following: JF(0) = const and 
asymptotically 

^ = W + o^ + -' xV>1, 

where a = 1/6. 

Again let us represent the integral 

7 ^ y V) tldtl 
J ^^ M ^( M 2 + 1) 5/2' 

as the sum of U\ and C/ 2 so that 



oo oo 



■u 3 



Now two first subtractions are needed to estimate U\ = X + Y + Z , where 

1 ^ C 1 ^ 

z = — y ^(m 2 ) Md« , y = - — — y j^(m 2 ) M 3 d« , 

x x 



o 



Clearly, Z is convergent at large x and approximately 



z « ±]r(u 2 )udu- JL. 

Further, a single counter-term is required for Y: 

_ 1 oo 

2 x 4 M 1 ; « + iy 2 x 4 

Now let us consider X. Here the series in the brackets begins with the term propor- 
tional to (w/x) 4 ) an d a single regulator is needed. The next term will require two and so 
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on. The total contribution of the order x 4 is given once JF is replaced by a/u A in X. The 
net result reads: 

a /37 5 , \ 



f/ 2 + x 



X 4 V12 
and, therefore, 

If,. , a f 31 5 , 1 N 

^ 2 = -/^>d. + ^- + -l g- 



5 " 



2x 4 •/ V " v ' u + 1 
o 

As the final result, the non-zero part of j is: 



du . 



P 1 (, A 5 1 i B\ 



where 



oo 

A = 3 J F(u 2 )udu «-2x 10" 1 , (52) 
o 

log4£ = — - 12 / ( m 3 ^(m 2 ) - - 1 - 1 du. 
15 / V 6(u + l)J 

In this case there are two corrections, and the lowest is of the type (rot l)± \ I x rot / 1 found 
in [10] for the current perpendicular to / (one should be referred to the formula (53) in 
[10]). The coefficient A (52) has been estimated numerically in [23]. The next term is the 
new cubic one and it includes the logarithm of the London parameter. 



5.3 EXAMPLE 3 

In this case we shall take into account the whole (37) which would imply appearance of the 
quadratic corrections of the type | / x rot / \ (v — (l/4m)rot/)|| and | / x rot / |(<9i/ 2 + d 2 h) 
[10]. However, we put here xi — Xi f° r simplicity, and so the answer expected would 
demonstrate such corrections only in principle. To this end we shall investigate the third 
component j C orr,3 which is along /. We obtain from (42): 

3p 1 ff dudv . „ . . „ . . 

((« - v)T- + (u + v)r+) , (53) 



2tt £oX 2 ^ ^J(u 2 + 1) 5 (1 -v 2 ) 



T T = j exp(-t(^) 2 Q) 



x (J_^i e xp((i - tanht)(xM) 2 (l-^ 2 )Q 3 ) I dt , (54) 
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where Q 2 stands for 

Q~ 2 = 1 + u 2 =f 2uv = 1 - v 2 + (ut v) 2 , 

the domain n is given by {(u, v) : u G [0, oo[, v G [0, 1]}, and u = tan#, v = cos0. 

The function Ai is still rather complicated and so the present consideration becomes 
less elegant then the two previous. The estimations we are interested in will be obtained 
without providing the asymptotic integral formulas for the coefficients. Besides, we shall 
assume that not only x but logx also is large (logarithmic accuracy), and thus only the 
logarithmic terms will be kept in the third order. Let us proceed estimating $ (40) in 
general situation. By steepest descent we get: 

*(x 2 .Q)--}^ + 2^, (55) 



at Q = \A\ /ac F > 1 (i.e. either A /ac F or x must be > 1, and it is forbidden to tend 
A 2 /ac F to zero). In the opposite case Q < 1 

$(x 2 , Q) ~ -1 + tt 1/2 x 2 (^f) V2 • (56) 

First of all, let us consider the contribution to j (53) which is due to u G [1, oo[. Here 
the function T can be expanded by steepest descent because (xu) 2 Q 3> 1. This expansion 
will begin with the third order term const x (^oX 6 )" 1 which is not of interest for us. So, 
in what follows we shall take < u < 1 in j. 

Now let us consider the domain < u < 1/x- Approximately we put: 

o 1 Vx i 

^ to? I * / 5^7^ «« - + (« + ^) ■ ("J 

' 

In this case A 2 /ac F ~ (x«) 2 (u =F ^) 2 and the Eq.(56) should be used as far as A 2 /ac F 
can become zero, while the main contribution is due to the region where A 2 /ac F is strictly 
less than 1. So we obtain 

j^-i + tt 1 / 2 *" (1 -y 2) . (58) 

The part of j which is due to the first term in (58) looks as follows: 

-^I^'l udu tjT^W = -^^- (59) 

The second term in (58) does not contribute at v > 1/x because sign(w — v) — —1, and 
therefore 

3p 1 VX VX p 1 

7Tl= 7 — udu dv(sign(w - u) + 1) = —= — — . (60) 
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Then, the total contribution below u — 1/x is: 



' - hM£ - 4) ■ 

At last, let us consider the rectangle {(u, f) : 1/x < u < 1, t) € [0, 1]}. Here F + can 
safely be expanded by the Laplace method. As to JF_, the integral diverges when A 2 /ac F 
falls into the strip \u — v\ < 1/x along the diagonale u = v but the resulting singularity is 

integrable. Here there is no interesting contribution as far as x 2 > 1 and the strip's width 

is 2/x- Outside the strip the use of (55) allows to fix unambiguously the logarithmic third 
order term. The experience of the previous calculations shows us that the coefficient at 
(£oX 4 ) -1 is mainly due to < u < 1/x, and therefore its order of magnitude should be 
given by (61). So, we get with JF + and JF_: 



p 1 f du 1 f dv 

i/x 



f du 1 r dv / 2 2 \ 

y "^( M 2 + 1) 5/2 y (i_ V 2)i/2 [ u + 8u - 5 j • 



The last equation results in the following contribution: 

8 6x 6 X 

while the total quadratic correction (at least the order) is given by 

_3\ ^_ 37xl(r 2_ (63) 

2 £oX 4 V 2 y 2 £oX 4 1 ; 



6 DISCUSSION 

The present paper is concerned with the two main problems: to calculate the mass cur- 
rent j in weakly inhomogeneous He — A using thermal Green functions and to obtain 
its asymptotic expansions at T = provided the London limit holds. Two main assump- 
tions are of importance for our approach: the static order parameter can be linearized 
due to slowness of its spatial variation, and only those first order differentiations are 
retained in the chosen mixed representation which are due to the kinetic energy of the 
BCS— Hamiltonian. Using slowness of the orbital vector texture we reduce the three ini- 
tial dimensions to the one-dimensional situation so that the resulting operator 7i em has 
the simple form of the Hamiltonian of the Landau problem. Therefore, one can solve the 
governing Dyson-Gorkov equation exactly: just using the eigen-functions of H e m [22-24]. 
Thus a collection of exact formulas both for the Green function and j appears which opens 
the possibility to derive systematically the higher gradient corrections to the dominant 
expression j (1). The present paper completes the preceding ones [22-24] which have 
been aimed at a more thorough resolution of the mathematical difficulties of the Ref. [10]. 
Our approach unravels the situation and provides a correct procedure to find the struc- 
ture and the order of magnitude of higher contributions to (1). The given approach is 
manifestly advantageous because the Laplace method is highly appropriate in the London 
limit. 
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Mathematically, we are mainly concerned with the ^-integrated and then cu-summated 

— # 

normal Green function which results in two representations for j: the series and the 
integral. The integral one seems to be more attractive as far as it allows to obtain a self- 
contained expression (41) for the net correcting contribution. The last can satisfactory 
be studied by steepest descent due to the London limit holds, i.e. it can be expanded in 
powers of the orbital vector gradients. Particular limiting cases (zero temperature limit 
following the limit of lowest order in gradients, and vice versa) confirm the correctness of 
our manipulations. The known expression (1) arises as the lowest approximation without 
resorting to specific gauges for the order parameter. The appropriate "commutative" 
diagram is considered in the Sec. 4. 2. 

Three special orientations of rot / are considered in the Sec. 5 to deduce the correcting 
terms explicitly: rot / is parallel (Example 1), and perpendicular (Example 2) to /, while 
the third Example 3 implies an intermediate orientation of rot / with respect to /. Cor- 
rections are considered up to third order in gradients of / and they are not only of pure 
polynomial type. Namely, new cubic corrections are found which contain the logarithm of 
the London parameter. By comparison with [10], we provide the numerical coefficients at 

the second order terms. The nodes vicinities < 9 (or n — 6) < 1/x on the Fermi sphere 

give the dominant contribution to the numerical coefficients at the pure power terms. 

In the first case only the third order logarithmic correction is present. In the second 
case both the lowest ones appear: the quadratic and the cubic. As it is clear from the 
analysis [10], all the quadratic corrections should be proportional to | / x rot / 1, and that 
correlates with the absence of quadratic term in Example 1. We have specified the Example 
2 so that jcorr is orthogonal to /. Thus the second order term in (51) would correspond to 
that which has been written in [10] in the form (rot/)j_| / x rot / 1, and the corresponding 
numerical coefficient is A (52) [23] . Example 3 also results in the both lowest corrections: 
the quadratic (63) and the cubic (62). As to the quadratic correction along /, the answer 
of [10] reads: 

p£ \T x rotf| (A(v 3 - rot //4m) + B (dj 2 + d 2 h)) ■ 

The term (63) above is just to be compared with the last expression at A ^ 0, B = 0. 
It is clear from (36) that the result obtained by us is also applicable to establish the 
contribution at A = 0, B ^ 0. 

Moreover, our situation is more rich than in [10] because the logarithmic corrections 
are demonstrated. Indeed, let us recall the correction found in [13]. As it is seen from 
(36), the components of the superfluid velocity mv and gradients of / enter equally to 
the parameter a, and therefore Xorb (v s • d) l a d l a , where Xorb is logarithmically large [9], 
should be treated as the logarithmic third order term, and it would correspond to the 
third order contribution in (62). 

To conclude, the investigation presented which is based on thermal Green functions 
demonstrates various corrections of second and third order to the mass current (1) at 
T = provided the London limit condition holds. The corrections predicted early in 
[10, 13] can systematically be deduced in our approach. The representations obtained for 
j and j corr would serve a basis for further investigations. 
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APPENDIX I 

The following Mehler formula is valid: 

00 1 / 1 — n \ 

^0 ^(l-a 2 ) V l + a ) 

where ^ n (y) are the Chebyshev-Hermite functions (see generating functions of Hermite 
polynomials in [28]). Using 

00 

(n + q)- 1 = J dte- t(n+q) , n > 0, q > 0, 


and (All), we obtain: 

S^-v!r/OT-((T-«)—'^) 

[22]. 

APPENDIX II 

Here we obtain a (35) as in [10]. To this end let us rewrite 

A(k, r) = 5{k ■ Ai(r) + ik ■ A 2 (f)) (AII.l) 

with the help of 

A Q (f) » A (O) + r ■ dA a (0) , (a = 1, 2), (AII.2) 
where the derivatives diA a (0) are linear functions of A a (0) : 

d t A a (0) = u? x A a (0). (AII.3) 

The Eq. (AII.3) acquires the form of the Mermin-Ho relation [30] provided the identifica- 
tions 

= -2mvi , c4 = dih , u\ = -dit 2 (AII.4) 
18 



are made. Substituting (AII.2)-(AII.4) to (AII.l) one gets 

A(k,r) = 5(%i + ik 2 ) + piS (2mv p k p (ki + i% 2 ) — k 3 k p (d p l2 — id p l] 

= 5 sinfle** + p [<*(...) e^' 2 '^] (AII.5) 

Eventually, the square brackets in (AII.5) are denoted as a, and the factor e 1 ^ is to make 
it a real positive fixing therefore A\ , A 2 in the plane perpendicular to I. From (AII.5) it 
is seen that acF is a linear form of gradients which can be written formally as 

.^^(pdWl^S^. (A...6) 
By (8) we consider acp/S 2 as small parameter. 
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